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ABSTRACT 

We present absolutely calibrated far-ultraviolet spectrophotometry of the quasar 3C 273 covering the 900-1800 
A range. Our ~ 3 A resolution spectra were obtained with the Hopkins Ultraviolet Telescope during the Astro- 
1 mission in December 1990 and during the Astro-2 mission in March 1995. Both spectra exhibit a change in 
slope near the Lyman limit in the quasar rest frame. At longer UV wavelengths, the continuum has a power-law 
index of 0.5-0.7, while shortward of the Lyman limit it is 1.2—1.7. The energy distribution in vf v therefore peaks 
close to the quasar Lyman limit. The short wavelength UV power law extrapolates well to match the soft X-ray 
excess seen in simultaneous observations with the Broad-Band X-ray Telescope and nearly simultaneous ROSAT 
observations. 

The general shape of the broad-band spectrum of 3C 273 is consistent with that of an optically thick accretion 
disk whose emergent spectrum has been Comptonized by a hot medium. Our UV spectrum is well described by 
a Schwarzschild black hole of 7 x 10 8 M Q accreting matter at a rate of 13 M Q yr" 1 through a disk inclined at 60 
degrees. Superposed on the intrinsic disk spectrum is an empirically determined Lyman edge of optical depth 0.5. 
The Comptonizing medium has a Compton parameter y w 1, obtained with an optical depth to electron scattering 
of unity and a temperature of 4 x 10 s K. 

This overall shape is the same as that found by Zheng et al. and Laor et al. in their UV and X-ray composite 
spectra for quasars, giving physical validity to the composite spectrum approach. When combined with those 
results, we find that the generic ionizing continuum shape for quasars is a power law of energy index 1.7-2.2, 
extending from the Lyman limit to ~1 keV. The observational gap in the extreme ultraviolet for these combined 
data describing the quasar continuum shape is now only half a decade in frequency. 

Subject headings: Accretion Disks — Galaxies: Quasars: General — Galaxies: Quasars: Individual (3C 273) — 
Ultraviolet: Galaxies 



1. INTRODUCTION 

It has long been known that the energy distribution of 
quasars, expressed as vL u , must have a peak somewhere in the 
far-ultraviolet to soft X-ray region of the spectrum. Shields 
(1978) first suggested that the optical and near UV flux of 
quasars might be due to the Rayleigh-Jeans portion of a black 
body spectrum whose peak lies in the unobserved extreme UV 
region, arising from an op tically thick accreti on disk surround- 
ing a massive black hole ( Lynden-Bell 1969 ), and dubbed the 
"big blue bump." This idea was then further developed and ap- 
plied to UV obs ervations of quasars obtained with IUE (M alkan 
& Sargent 1982; Malkan 1983| ). It was subsequently found that 
the X-ray spectra of quasars also showed a "soft excess" of flux, 
above the extrapolation of the relativel y flat power-law that fit 



and the soft X-ray (10 17 Hz), and in particular, establishing 
more precisely wher e in the extreme UV the peak of the big 



blue bump is located (|Arnaud et al. 19851 ; [Bechtold et al. 1987 
O'Brien, Gondhalekar, & Wilson 19881). 



the spectrum at ene rgies above 2 keV (Turner & Pounds 1988 



Masnou et al. 1992), and which might arise from the Wien por- 



tion of the same thermal spectrum. In addition to whatever light 
this spectral region might shed on the black-hole accretion-disk 
model for quasars, it is also a fundamental input to studies of 
t he photoionization of the broad emission line clouds in quasars 
( Krolik & Kallman 1988) as well as the photo ionization of the 
intergalactic medium ( Haardt & Madau 1996| ). Consequently 
there has been widespread interest in determining the detailed 
shape of the spectrum of quasars between the optical (10 15 Hz) 



The quasar 3C 273 provides an excellent candidate for ad- 
dressing this problem because of its low redshift (z=0.158) and 
its high flux at both UV and X-ray wavelengths. I t was the first 
quasar detected in X-rays ( Bowyer et al. 1970 ) and also the 
first quasar observed in the far-ultraviolet (D avidsen, Hartig, & 
Fastie 1977), both detections having been achieved with sound- 
ing rocket experiments before the advent of effective satellite 
observatories. 

At X-ray energies the spectrum of 3C 273 in the 2-10 keV 
range is well fit by a power-la w, F u oc v~ a , with energy index 
a w 0.5 ( Worrall et al. 1979 ), bu t observations at low er en- 
ergies with EXOS AT and Ginga (Turner et al. 1985[ C our- 
voisier et al. 1987, [Turner et al. 1990[), Einstein (Wiik es & 
Elvis 19 87, [Turner et al. 199U Masnou et al. 1992[ ), and 
ROS AT (IStaubert et al. 1992|; L each, McHardy, & Papadakis 
1995, hereafter LMP95) have established the existence of a soft 
excess. Reporting on simultaneous observations of 3C 273 in 
1990 December with ROSAT (0.1-2 keV) and Ginga (2-10 
keV), Staubert et al. (1992) obtained a good fit with the sum of 
two power laws, with a = 0.56 dominating at high energies and 
a = 2.5 at low energies. However, the soft component could 
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Table 1 
Observation Log 





Instrument 


Date 
Start Time (GMT) 


Exposure 
(sec) 


Comment 


Astro- 1 


HUT 


1990 December 05 

23:28:24 
1990 December 09 

17:23:54 


1251 
2138 






BBXRT 


1990 December 09 
17:37:01 


1598 






IUESWP40331 


1990 December 14 


1800 






SWP40333 




1500 






T WP1 941 1 




1800 






LWP19412 




1800 






SWP40391 


1990 December 19 


1800 






SWP40392 




1800 






LWP19447 




1800 






LWP 19448 




1620 






ROSAT 


1990 December 18-21 


497 






KPN0 2.ini 


1990 December 05 


60 




AA 










Astro-2 


HUT 


1995 March 12 
06:05:03 


1180 


Flux correction 1.02 






1995 March 15 


1304 


Flux correction 1.05 






07:11:43 







also be modeled with a black body s pectrum or with thermal 
bremsstrahlung ( jStaubert et al. 1992 ). More extensive obser- 
vations of 3C 273 with ROSAT have been reported by LMP95. 
They established that "two physically distinct emission com- 
ponents are present", with the spectrum best modeled by a 
combination of two power-laws with absorption by the Galac- 
tic column of gas, Nh = 1-84 x 10 20 cm" 2 . Constraining the 
harder component to have a/, = 0.5 as observed by EXOSAT 
and Ginga, the soft component was found to be significantly 
steeper, with a s = 1.7 for most of the observations. Attempts 
to fit the soft component with various other models, including 
blackbody and bremsstrahlung radiation were found unaccept- 
able (LMP95). 

At ultraviolet wavel engths the first determination of the spec- 
tral index of 3C 273 (Davidsen, Hartig, & Fastie 1977) gave 
a = 0.6 over the range from the optical to Lyman a, with a broad 
bump around 3500 A, perhaps associated wit h Balmer contin- 
uum and Fe II line emission (Baldwin 1975), that has some- 



times been called the "small blue bump". It was also pointed 
out by Davidsen et al. that the spectral index must increase 
to a > 1.2 somewhere near or above the Lyman limit in or- 
der to agree with the X-ray observations. Subsequent exten- 
sive observations with IUE give a = 0.66 between 3100 and 
1250 A ( [Zheng & Malkan 1993] ). A very low resolution ob- 
servation with the Voyager UV Spectrometer in the 900-1200 
A band suggested the possibility of a spectral break at about 
1000 A (Reichert et al. 1988), perhaps associated with the Ly- 
man edge. Recent observations in the 900-1200 A region with 
ORFEUS imply a turnover in the spectrum of 3C 273 near the 
IUE/ORFEUS tr ansition wavelength at about 1200 A ( Appen- 
zelleretal. 1998). 

Here we report the results of observations of 3C 273 in the 
900-1800 A region made with the Hopkins Ultraviolet Tele- 
scope (HUT) on both the Astro- 1 and Astro-2 space shuttle mis- 
sions. The Astro- 1 far-ultraviolet data are accompanied by con- 



temporaneous observations at optical, ultraviolet, and soft and 
hard X-rays, yielding a quasi-simultaneous broad-band spec- 
trum for this quasar extending over more than 3 decades of fre- 
quency. The HUT data alone show a definite spectral break, 
with a peak in the energy spectrum vL v at about 920 A in the 
rest frame. The best-fit power-law at shorter wavelengths from 
920 A to 787 A, when extrapolated, nearly matches the simulta- 
neous soft X-ray spectrum. The Astro-2 HUT data are of even 
higher quality than those obtained on Astro- 1 and yield sim- 
ilar results, confirming that the peak of the "big blue bump" 
has now definitely been seen in 3C 273. A comparison of these 
data with a composite quasar spectrum in the ultraviolet (Z heng 
et al. 1997) and the X-ray bands ( |Laor et al. 1997 ), suggests 
strongly that a peak in the energy spectrum of all quasars occurs 
near the Lyman limit, and that the ionizing continuum from the 
Lyman limit to 1 keV is a power-law of a = 1 .7 - 2.2. 

2. OBSERVATIONS 

The Hopkins Ultraviolet Telescope (HUT) incorporates a 
0.9-m primary mirror and a prime-focus Rowland circle spec- 
trograph as described by Davidsen et al. (1992) and Rruk et 
al. (1995). First-order spectra recorded by the photon counting 
detector cover the 820-1840 A spectral range with a sampling 
of about 0.52 A per pixel and a point source resolution of about 
3 A. Absorption by interstellar hydrogen limits the observed 
spectra to wavelengths longer than the Lyman limit at 912 A 
(787 A in the quasar rest frame). 

HUT was used to observe 3C 273 on the Astro- 1 space shut- 
tle mission in 1990 December and again on the Astro-2 mission 
in 1995 March. The performance and calibration of the instru- 
ment are described for Astro- 1 by Davidsen et al. (1992) and 
Kruk et al. (1997) and for Astro-2 by Kruk et al. (1995, 1999). 
Because substantial changes to HUT's performance were made 
between the two flights, the results of the spectrophotometric 
observations may be regarded as independent measurements, 
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FIG. 1 . — 3C 273 spectrum (histogram) from 900 to 1820 A as observed with HUT on Astro-1 (lower panel) and Astro-2 (upper panel). Pixels are about 0.51 A 
and the resolution is about 3 A. Curves at the bottom of each panel give the values of the error array produced by the standard HUT reduction process. Strong peaks 
in the ltr error array for Astro-1 correspond to strong geocoronal emission lines (airglow) for this observation, obtained mostly in daylight. Airglow lines have been 
subtracted using templates developed from blank sky observations, and their positions are marked with ffi. The Astro-2 spectrum was obtained at night, resulting in 
much weaker airglow. Emission and absorption lines are marked. 



effectively produced by two different instruments. The Astro- 
1 observations were (necessarily) conducted primarily during 
daylight portions of the shuttle orbit, where there is substan- 
tial contamination by geocoronal emission lines that must be 
subtracted in the data reduction process. The Astro-2 obser- 
vations, however, were obtained mostly during night-time por- 
tions of the orbit, where the geocoronal contamination is greatly 
reduced. Only the night portion of the Astro-2 data is used. 
Details of the observations are given as an observation log in 
Table 1. 

For Astro-1 there are several nearly simultaneous observa- 
tions that can be used to determine the broad-band spectrum of 
3C 273. Also mounted on the shuttle with HUT was the Broad 



Band X-Ray Telescope (BBXRT) ( |Serlemitsos et al. 1992[ ), 
which observed 3C 273 in the 0.3-12 keV band with a res- 
olution of about 100 eV. ROSAT observations in the 0.1-2.4 
keV band were obtained less than two weeks after the Astro-1 



observations (Staubert et al. 1992), and Ginga observations in 
the 2 -20 keV band w ere also obtained during 1990 December 



( [Staubert et al. 19"92 ). In the ultraviolet, IUE spectra in both 
the short and long wavelength cameras were obtained between 
one and two weeks after the HUT observations, and these have 
been extracted from the IUE archive. Finally, an optical spec- 
trum of 3C 273 was kindly obtained for us by R. Green on the 
KPNO 2.1 m telescope during the Astro-1 mission. All these 



contemporaneous measurements are listed in the observation 
log in Table 1 . 

The HUT observations of 3C 273 from Astro-2 in 1995 
March are also listed in Table 1 . Unfortunately, we have been 
unable to find any nearly simultaneous observations at other 
wavelengths for this epoch. 

3. RESULTS FROM THE HUT OBSERVATIONS 
3.1. Observed Spectra 

The two HUT observations of 3C 273 from Astro-1 were 
summed and then processed and reduced to an absolute flux 
scale following our standard procedures as described by Kruk 
et al. (1997). For Astro-2 the two observations were reduced 
separately because of th e variation of the i nstrument sensitiv- 
ity during the mission (Kruk et al. 1995), and the resulting 
spectra were then averaged, weighted by their exposure time, 
to yield the final spectrum presented here. As noted in Table 1, 
the Astro-2 data had slight flux corrections of 2-5% applied due 
to light loss at the slit induced by pointing errors. The results 
for both missions are displayed in Fig. 1. 

It is worth emphasizing that the Astro- 1 and Astro-2 spectra 
were obtained, reduced, and calibrated under markedly differ- 
ent circumstances. The Astro-1 data were obtained with an in- 
strument of relatively modest sensitivity, almost entirely during 
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orbital day when there is substantial contamination by geocoro- 
nal emission lines that must be subtracted, and were absolutely 
calibrated by reference to an observation of the DA white dwarf 
G191-B2B u sing a model atmosphe re computed for this star by 
P. Bergeron (Davidsen et al. 1992). The Astro-2 data, on the 
other hand, were obtained with an instrument of much higher 
sensitivity, entirely during orbital night with minimal geocoro- 
nal contamination, and were absolutely calibrated by reference 
primarily to observations of the DA white dwarf HZ43 using a 
model atmosphere compu ted by D. Finley u sing D. Koester's 
model atmosphere codes (Kruk et al. 1995). In spite of these 
major differences, the resulting spectra in Fig. 1 are remarkably 
similar. Of course, 3C 273 is known to vary, so the close agree- 
ment between observations made more than four years apart is 
fortuitous. Indeed there are clear changes observed both in the 
continuum and in the emission line strengths (discussed below). 
However, both spectra show a change of continuum slope near 
the Lyman limit in the rest frame of the quasar, which is further 
discussed in the next section. 

3.2. Empirical Models for the Spectra 

We first remove the airglow lines from the spectra. We use a 
geocoronal Lya template that is derived from blank fields and 
scaled by the exposure time, and subtract it from the HUT spec- 
trum of 3C 273. Other airglow lines are fitted and removed, 
but significant residuals still exist in the regions around the 
strongest airglow lines at 1216, 13 02 and 989 A . 

We use the IRAF task specfit ((Kriss 1994]) to fit both the 



Astro- 1 and Astro-2 spectra with various components for the 
continuum, emission lines and absorption features. We use the 
same overall model for both spectra, and identical fitting win- 
dows spanning the wavelength intervals 912-1 194, 1238-1287, 
and 1319-1820 A. For the Astro-1 data, we also model the re- 
maining airglow residuals near 989 A. 

The intrinsic emission lines of Lya, C IV and O VI are mod- 
eled with dual Gaussians, and other emission components are 
modeled with single Gaussians. We constrain many compo- 
nents of the fit by tying related parameters together. The full- 
widths at half maximum (FWHM) of the weaker broad emis- 
sion lines (S VI, C III, N III, O I, C II) are all linked, and their 
wavelengths are linked to that of narrow Lya by the ratio of 
their laboratory wavelengths. The FWHM of broad O VI A1034 
is linked to broad C IV A1549. The wavelengths and FWHM of 
the N V components are linked to the corresponding C IV com- 
ponents. We include one unidentified weak, broad emission 
feature at a rest wavelength of 1074 A. This is not He II A1085, 
and it is approximately at the location of the unidentified fea- 
ture also seen in Faint Object Spectrograph (FOS) observations 
byLaoretal. (1995). 

For the absorption lines, we start by including features iden- 
tified in pr evious, higher resolution observations (FOS: B ahcall 
etal. 19 91; GHRS: Morris et al. 199U ORFEUS: H urwitz et 
al. 1998). We then add additional features as required to ob- 
tain a good fit to the HUT data. As with the emission lines, we 
constrain many components of the fit by tying the FWHM of 
the lines together in groups by wavelength regions where the 
instrumental resolution is roughly the same. All lines are as- 
sumed to be unresolved, except for the blend of O VI, C II, and 
O I around 1038 A, the blended Si II lines at 1192 A, and the 
blended C IV doublet at 1549 A. 

The wavelength region below 950 A is heavily absorbed 
by the high-order lines of Galactic hydrogen. To model this, 



we use a set of absorption profiles, each including 50 Lyman- 
series lines, with a fixed column density of 1.8 x 10 20 cm -2 . 
The absorption lines are modeled using Voigt profiles cover- 
ing a range of 10-20 kms -1 in the Doppler parameter which 
are then convolved with the instrument resolution, a Gaus- 
sian of 3 A FWHM. For our extinction correction, we adopt 
E(B-V) = 0.032, as established by Lockman & Savage (1995), 
and use the extinction curve of Cardelli, Clayton & Mathis 
(1989), with Ry = 3.1. This value of the extinction is also con- 
sistent with the E(B-V) = 0.03 adopted by Lichti et al. (1995) 
and Appenzeller et al. (1998). 

For both the Astro- 1 and Astro-2 spectra, we obtain our best 
fit using a broken power-law model for the continuum shape. To 
ascertain the significance of the apparent spectral break near the 
Lyman limit, we compare fits performed using a single power- 
law for the continuum shape to a broken power-law. The sin- 
gle power-law fits are summarized in Table 2, and the broken 
power-law fits are summarized in Table 3. For the Astro-1 data, 
the fits with a single power law do not produce satisfactory re- 
sults. These fits are shown as smooth curves in Fig. 2 and Fig. 3. 
Note that the single power-law fit deviates significantly at wave- 
lengths shortward of 1200 A. Even allowing for reddening cor- 
rections higher than the Galactic value, the fits are still poor, as 
shown in Table 2. The spectral break in the sub-Lya region is 
too sharp to be accounted for by reddening. 



TABLE 2 
Single Power-Law Fits to 3C 273 



E(B-V) 


Astro- 1 




Astro-2 




a 


X 2 /dof 


a 


X 2 /dof 


0.01 


0.94 


1689/1544 


1.16 


1723/1501 


0.02 


0.81 


1651/1544 


1.00 


1698/1501 


0.032 


0.62 


1619/1544 


0.79 


1697/1501 


0.04 


0.49 


1615/1544 


0.67 


1708/1501 


0.05 


0.34 


1602/1544 


0.50 


1746/1501 


0.06 


0.18 


1607/1544 






0.08 


-0.15 


1650/1544 







TABLE 3 

Broken Power-Law Fits to 3C 273 a 



Parameter 


Astro-1 


Astro-2 


,b 
Jo 


4.37 ±0.04 


4.41 ±0.02 


Ao (A) 


1064 ±8 


1036 ±6 


ai 


1.69±0.16 


1.19±0.15 


OC2 


0.46 ±0.03 


0.74 ±0.02 


X 2 /dof 


1584/1542 


1692/1499 



"Assumes extinction fixed at E(B — V) = 0.032. 

b Measured at the break wavelength Ao, in units of 10* 13 



ergs s 1 cm 2 . 



Although the Astro-2 spectrum is fairly similar to the Astro-1 
data, the continuum shape is noticeably less blue at wavelengths 
longward of Lya, and the broad emission lines are significantly 
brighter. One can see this most easily in the comparison of the 
fitted spectra shown in Fig. 4. Since the Astro-2 spectrum is not 
as blue, the break at shorter wavelengths is not as dramatic, and 
a single power law reddened at the nominal Galactic value fits 
only slightly worse than the broken power law as one can see 
in the top panels of Figures 2 and 3. Thus, we conclude that a 
short wavelength spectral break is required by our Astro-1 data, 
but not by the Astro-2 data. 
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FIG. 2. — Observed spectra of 3C 273 with the best-fitting empirical models 
as described in Section 3.2. The models include emission and absorption lines 
and a continuum described either by a single power-law (dotted line) or by a 
broken power-law (solid lines). Reddening with E(B — V) = 0.032 is applied to 
the models. Broken power-laws provide the best fit for both spectra. A single 
power-law gives an unacceptable fit for the Astro- 1 spectrum; for Astro-2 a sin- 
gle power-law is worse, but acceptable. Parameters of the models are given in 
Tables 2-7. 



error matrix of the fit, and they represen t the 1 a con fidence in- 
terval for a single interesting parameter (Kriss 1994). 
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FIG. 3. — Same as Fig. 2 but showing the crucial 900-1200 A region at an en- 
larged scale. Positions of the Lyman series of interstellar hydrogen absorption 
lines are indicated. Series lines up to L7 are clearly seen in the data. The con- 
tinuum in this region is well-fit by the broken power-law model (solid curves), 
while the best single power-law model (dotted) lies above the data at the shortest 
wavelengths. 



As the broken power law continuum models give our best fits 
for both data sets, we quote the parameters for the emission and 
absorption lines using these continuum models. Tables 4 and 
5 give the observed wavelengths, the observed fluxes, and full- 
widths at half-maximum (FWHM) for the Astro- 1 and Astro-2 
observations of the emission lines, respectively. Tables 6 and 
7 list the observed wavelengths, equivalent widths (W\), and 
FWHM for the absorption lines. None of the line widths are 
corrected for the instrumental resolution. For all absorption 
lines, the best-fit widths are consistent with the instrumental 
resolution and should therefore be considered unresolved. The 
error bars quoted on all parameters are determined from the 



4. THE BROAD-BAND SPECTRUM OF 3C 273 

The BBXRT data were obtained simultaneously with the 
Astro- 1 HUT data, and we have retrieved them from the Na- 
tional Space Science Data Center (NSSDC) archive. To obtain 
the unfold ed X-ray spec trum, we use the data reduction pack- 
age xspec ( Arnaud 1996j ). The data points in channels 1-20 and 
450-512 are excluded because of high uncertainties. We bin 
the data so that each bin has at least 20 counts. Our preliminary 
fits with single power laws or broken power laws yielded poor 
fits, with a significant deviation around 0.6 keV. Done (1993) 
suggests that this is a blueshifted O VIII resonance absorption 
feature. Empirically modeling this feature as an absorption 
edge, we use a dual power-law with fixed Galactic absorption 
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to model the X-ray spectrum: 



dN/dE = +f 2 E-^)exp(-T(E/E, h y 3 )exp(-T(N H ,E)). 

This model gives a good fit (% 2 = 52.5 for 52 binned data 
points and 6 free parameters). The parameters for our best 
fit are: the photon power-law index a\ = 4.3 ± 0.9, a 2 = 
1.64 ±0.09; flux at 1 keV: f x = 0.0081 ± 0.0030 and f 2 = 
0.017±0. 0025 photon s" 1 cm" 2 keV _1 ; and the absorption edge 
at E th = 0.5 1 ± 0.03 keV with an optical depth r = 1 . 17 ± 0.3 1 . 
This fit and the residuals are shown in Fig. 5. 
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FIG. 4. — Direct comparison of the fits to the HUT Astro- 1 (dotted lines) and 
Astro-2 (solid lines) spectra. Although the continuum fluxes are roughly the 
same, there are clear differences in the line fluxes and continuum shapes. 

The nearly simultaneous optical, ultraviolet, and X-ray data 
collected in 1990 December during or shortly after the Astro- 1 
mission were listed in Table 1 and are shown in Fig. 6. The 
BBXRT and ROSAT data shown in Fig. 6 are corrected for 
Galactic absorption, using a Morrison & McCammon (1983) 
model with N H = 1 -8 x 10 20 cm" 2 . (Note that we show unfolded 
X-ray data in Fig. 6, and that the plotted fluxes are model de- 
pendent.) Also shown is the power-law that best fits the short 
wavelength portion of the HUT spectrum as part of the broken 
power-law fit described in Section 3.2. This component has 
a = 1.7. Two additional power-laws are also shown as dashed 
lines with indices of 1.7 ±0.36, representing the 1-er errors 
from our fit added linearly to the uncertainty of ±0.2 corre- 
sponding to a range in E(B-V) of ±0.01, which is our estimate 
of the uncertainty in the extinction for 3C 273. It is clear from 
Fig. 6 and Fig. 2 that the energy distribution of 3C 273 peaks 
within the HUT spectral range, a little longward of the Lyman 
limit in the quasar rest frame, and that the short wavelength 
power-law that fits the HUT data extrapolates well to fit the soft 
X-ray data from ROSAT and BBXRT. 



There still remains a considerable gap of one decade in fre- 
quency that cannot be filled for 3C 273, however, because the 
ultraviolet spectrum cannot be extended past the Galactic Ly- 
man limit at 787 A in the quasar rest frame, and the soft X-ray 
spectrum is severely attenuated by interstellar absorption be- 
low 0.2 keV. Indeed this extreme ultraviolet gap cannot readily 
be decreased in the spectrum of any single quasar, since the 
higher redshifts that would allow the ultraviolet spectrum to be 
extended will at the same time limit the observed X-ray spec- 
trum to higher rest-energies. We therefore turn to a composite 
quasar spectrum assembled from many different objects to try 
to reduce the size of the extreme ultraviolet gap. 

In a previous paper (Zheng et al. 1997) we have presented the 
results of an analysis of HST FOS spectra of 101 quasars, most 
with redshifts from 0.33 to 1.5, with a few high redshift objects 
with z up to 3.6. The composite quasar spectrum was also fit 
with a broken power-law, the break occurring around 1050 A in 
the rest frame. The short wavelength component was found to 
have a = 2.0 for the total sample (which is subject to unknown 
selection effects), while the subset of 60 radio-loud quasars had 
a = 2.2 and the subset of 41 radio-quiet quasars had a = 1.8. 

In the lower panel of Fig. 6 we plot the radio-loud composite 
quasar spectrum from Zheng et al. (1997) scaled to match the 
flux observed for 3C 273 at 1000 A, near the break point in both 
spectra. Also shown as a solid line (with a break at 2 keV) is a 
composite quasar X-ray spectrum for radio-loud quasars from 
Laor et al. (1997), which has been scaled in flux to match our 
UV-optical composite using the average a ox of 1 .445 from the 
sample of Laor et al. 
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FIG. 5. — Upper panel: the best-fit model to the Astro-1 BBXRT X-ray spec- 
tram. As described in §4, the solid line is a dual power law model with Galac- 
tic absorption and an empirical absorption edge at 0.5 keV folded through the 
BBXRT response function. The error bars are \a. Lower panel: residuals to 
the fit. 
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TABLE 4 










Astro- 1 Emission Lines 






Feature 


Aote 


Flux" 


FWHM 




(A) 




(km s 


') 


S VI, A 933.38 


1079.42 ±0.20 


2.0±1.0 


4905 ± 


721 


S VI, A 944.52 


1092.32 ±0.20 


1.0±0.5 


4905 ± 


721 


C III, A 977.03 


1129.95 ±0.20 


3.8± 1.2 


4905 ± 


721 


N III, A 991.00 


1146. 12 ±0.20 


2.9 ±1.1 


4905 ± 


721 


O VI, A1034.00 


1191.46±0.20 


16.1 ±3.9 


4225 ± 


776 


(J VI, A1034.00 


1191.46±0.20 


13.7 ±5.5 


10554 ± 


723 


VI, A1034.00 b 




29.8 ±6.8 






No ID 


1243.66 ± 1.15 


5.3 ±1.1 


4905 ± 


721 


f ■ TTT \ 1 nr HC\ 

C ill, All ID. 10 


1360.28 ±0.20 


6.6± 1.8 


4905 ± 


721 


Lya, A 1213.0/ 


1406.21 ±0.20 


56.3 ±3.8 


2782 ± 


113 


Lya, Alzo.o/ 


1398.72±1.66 


49.8±5.1 


8375 ± 


733 


Lya, A1215.67 b 




106.1 ±6.3 






N V, A1240.15 


OA 1 ACT 

1431.80±0.32 


6.1 ± 1.8 


2782 ± 


308 


N V, A1240.15 


1434.68 ±1.89 


41.0±3.8 


10554 ± 


723 


NT V \ I OAO 1 
IN V, AlZ'+U. ID 




47.1 ±4.2 






I, A1304.35 


1508.05 ±0.20 


5.8 ±1.1 


4905 ± 


721 


CII, A1335.30 


1544.34 ±0.20 


2.1 ±1.2 


4905 ± 


721 


Si IV, A 1402.77 


1620.89 ±1.62 


15.1 ±2.4 


6978 ± 1325 


C IV, A1549 


1789.44 ±0.52 


17.8±2.5 


2782 ± 


308 


C IV, A 1549 


1793.03 ±1.89 


36.1 ±4.1 


10554 ± 


723 


C IV, A1549 b 




53.9±4.8 







a l(r 13 ergcnr 2 s~> 

b Total flux for the narrow and broad components. 



TABLE 6 



ASTRO- 


1 Absorption Lines 






Feature 




W X 


FWHM" 




(A) 


(A) 


(km s 


') 


O I, A 948.69 


947.09 ±0.41 


0.4 ±0.2 


820 ± 


162 


N I blend, A 953.00 


953.28 ±0.41 


0.4 ±0.2 


820 ± 


162 


P II, N I, A 964.00 


964.94 ±0.41 


0.5±0.1 


820 ± 


162 


C III, A 977.03 


975.68 ±0.30 


0.7±0.1 


820 ± 


162 


S III, A1012.50 




< 0.2 b 


1147 


Si II, A1020.70 


1020.70 ±0.30 


0.1 ±0.1 


1162± 


128 


OVI, A1031.93 


1031.93 ±0.29 


0.6 ±0.1 


1162± 


128 


O VI, C II, O I, A1038.00 


1037. 19 ±0.29 


1.5 ±0.1 


1344 ± 


148 


Ar I, A1048.22 


1049.04 ±0.39 


0.5 ±0.1 


1162± 


128 


Fell, A1063.18 


1063.18±0.39 


0.5±0.1 


1162± 


128 


Ar I, A 1066.66 


1067.51 ±0.39 


0.5±0.1 


1162± 


128 


Nil, A 1084. 19 


1084. 19 ±0.72 


0.5 ±0.2 


1162± 


128 


Fell, A1096.88 




<0.2 b 


779 




Fell, Al 121.97 


1121.97±0.72 


0.3±0.1 


779 ± 


127 


NI, Al 134.63 


1134.96 ±0.72 


0.4 ±0.2 


779 ± 


127 


Fell, Al 143.22 


1142.81 ±0.72 


0.4±0.1 


779 ± 


127 


Si II blend, A 1192.00 


1191.56±0.72 


1.1 ±0.3 


1536 ± 


497 


NV, A1238.82 




().2 C 


536 




NV, A1242.80 




().l c 


536 




S II, A1250.50 


1250.40 ±0.29 


0.5 ±0.1 


536 ± 


137 


S II, Si II, A1260.00 


1260.43 ±0.67 


0.4 ±0.1 


536 ± 


137 


HI, A1275.19 




0.1 c 


536 




C II, A1335.30 


1335.19±0.34 


1.0 ±0.2 


964 ± 


176 


HI, A1361.42 




().l c 


964 




Ni II, A1370.09 




0.1 c 


589 




Si IV, A1393.76 


1393.35 ±0.51 


0.2±0.1 


589 ± 


100 


Si IV, A1402.77 


1402.31 ±0.51 


0.1 ±0.1 


589 ± 


100 


Sill, A1527.17 


1526.25 ±0.56 


0.7 ±0.2 


826 ± 


225 


CIV, A1549.50 


1548.73 ±0.26 


1.3 ±0.2 


715± 


137 


Al II, A1670.79 




<0.4 b 


554 





"All tabulated widths are consistent with the instrumental resolution, and all 
features should be considered unresolved. 

b The upper limit for this line assumed a feature at the nominal wavelength with 
a fixed FWHM as shown. 

c Our fits included a feature at the nominal wavelength with the EW and FWHM 
fixed at the values shown. 



TABLE 5 



Astro-2 Emission Lines 


Feature 


Aofw 


Flux" 


FWHM 




(A) 




(km s 


') 


S VI, A 933.38 


1079.59 ±0.15 


4.2 ±0.2 


5145 ± 


151 


S VI, A 944.52 


1092.47 ±0.15 


2.1 ±0.1 


5145 ± 


151 


C III, A 977.03 


1130.05±0.15 


7.4 ±0.5 


5145 ± 


151 


N III, A 991.00 


1146.20 ±0.15 


6.5±0.8 


5145 ± 


151 


OVI, A1034.00 


1191.12 ±1.18 


13.6±9.2 


2332 ± 


590 


OVI, A1034.00 


1191.12 ±1.18 


38.7±3.0 


10416 ± 


477 


OVI, A1034.00 b 




52.3 ±9.6 






No ID 


1243. 66± 1.15 


5.6±0.8 


5145 ± 


151 


CIII, Al 175.70 


1360.08 ±0.15 


10.0± 1.0 


5145 ± 


151 


Lya, A1215.67 


1405.94 ±0.15 


70.7 ±3.3 


2799 ± 


76 


Lya, A1215.67 


1399.49 ±0.70 


93.5 ±4.4 


8445 ± 


322 


Lya, A1215.67 b 




164.2±5.5 






N V, A1240.15 


1433.61 ±0.4U 


6.6 ±1.6 


3070 ± 


213 


N V, A1240.15 


1433.79 ± 1.22 


52.8±3.5 


10416 ± 


477 


N V, A1240.15 b 




59.4±3.8 






OI, A1304.35 


1507.65 ±0.15 


6.5 ±0.9 


5145 ± 


151 


CII, A1335.30 


1543.90 ±0.15 


5.3 ± 1.3 


5145± 


151 


Si IV, A 1402.77 


1622.38 ±0.95 


21.5±2.0 


6800 ± 


671 


CIV, A1549 


1791.30±0.40 


28.6 ±2.8 


3070 ± 


213 


CIV, A 1549 


1791.53 ±1.22 


57.8±3.2 


10416 ± 


477 


CIV, A1549 b 




86.4±4.2 







"lO" 13 ergcnr 2 s~> 

b Total flux for the narrow and broad components. 



TABLE 7 



Astro-2 Absorption Lines 


Feature 


^obs 


W X 


FWHM" 




(A) 


(A) 


(km s 


') 


O I, A 948.69 


947.61 ±0.47 


0.4 ±0.1 


894 ± 


130 


N I blend, A 953.00 


953.28 ±0.22 


0.8±0.1 


894 ± 


130 


P II, N I, A 964.00 


964.04 ±0.44 


0.3±0.1 


894 ± 


130 


C III, A 977.03 


976.48±0.15 


1.0±0.1 


894 ± 


130 


Si II, A 989.87 


989.31 ±0.15 


0.5±0.1 


894 ± 


130 


N III, A 991.00 


990.44 ±0.15 


0.3±0.1 


894 ± 


130 


S III, A1012.50 


1011.80±0.47 


0.6 ±0.1 


1202 ± 


129 


Sill, A1020.70 


1020.29 ±0.62 


0.4 ±0.1 


1202 ± 


129 


OVI, A1031.93 


1031.00±0.41 


0.6 ±0.1 


1202 ± 


129 


OVI, CII, OI, A1038.00 


1037.43 ±0.21 


1.6±0.1 


1386 ± 


149 


Arl, A1048.22 


1051.70±0.69 


0.5±0.1 


1410 ± 


318 


Fell, A1063.18 


1063.74 ±0.67 


0.4 ±0.1 


1410± 


318 


Arl, A1066.66 




< 0.2 b 


1410 




Nil, A1084.19 


1083.46 ±0.52 


0.7 ±0.2 


1410± 


318 


Fe II, A1096.88 




< 0.2 b 


767 




Fell, Al 121.97 


1121.89 ±0.46 


0.2±0.1 


767 ± 


126 


NI, A 1134.63 


1134.83 ±0.18 


0.7 ±0.1 


767 ± 


126 


Fell, Al 143.22 




< 0.2 b 


767 




Si II blend, A 1192.00 


1191.33±0.76 


1.9 ±1.8 


1654 ± 


513 


NV, A1238.82 


1238.91 ±0.77 


0.2 ±0.1 


588 ± 


143 


NV, A1242.80 


1242.60 ±0.77 


0.1±0.1 


588 ± 


143 


S II, A1250.50 


1250.63 ±0.64 


0.2±0.1 


588 ± 


143 


S II, Si II, A 1260.00 


1260.75 ±0.20 


0.5 ±0.1 


588 ± 


143 


HI, A1275.19 




0.1 c 


588 




C II, A1335.30 


1334.34 ±0.21 


0.7 ±0.1 


651 ± 


103 


HI, A1361.42 




0.1 c 


651 




Ni II, A 1370.09 




0.1 c 


599 




Si IV, A1393.76 


1393.68±0.17 


0.5 ±0.1 


599 ± 


100 


Si IV, A1402.77 


1402.63 ±0.17 


0.2 ±0.1 


599 ± 


100 


No ID, A 1492.00 


1492.28 ±0.54 


0.4 ±0.1 


599 ± 


100 


No ID, A1522.00 


1522.85 ±0.25 


0.2 ±0.1 


246 ± 


92 


Sill, A1527.17 


1526.68±0.17 


0.7 ±0.1 


453 ± 


96 


No ID, A1533.00 


1531.85 ±0.25 


0.5 ±0.1 


453 ± 


96 


C IV, A1549.50 


1549.23 ±0.33 


1.0 ±0.2 


884 ± 


147 


AlH, A 1670.79 


1670.24 ±0.42 


0.6 ±0.2 


539 ± 


170 



"All tabulated widths are consistent with the instrumental resolution, and all 
features should be considered unresolved. 

b The upper limit for this line assumed a feature at the nominal wavelength with 
a fixed FWHM as shown. 

c Our fits included a feature at the nominal wavelength with the EW and FWHM 
fixed at the values shown. 
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FIG. 6. — Upper panel: Quasi-simultaneous broad-band spectrum of 3C 273 
from optical to X-rays obtained in December 1990, corrected to the rest 
frame. The ultraviolet and optical data have been corrected for reddening with 
E(B- V) = 0.032, and the ROSAT and BBXRT data have been corrected for in- 
terstellar absorption with Nfj = 1 . 84 X 1 20 cnr 2 . To show the continuum shape 
more clearly, the UV and optical data have been smoothed with a 9-pixel boxcar 
filter. The best-fit power law for the short wavelength HUT data with a = 1 .7 is 
indicated, along with 1 a uncertainties for the slope. Lower panel: Composite 
radio-loud quasar spectrum in the UV-optical from Zheng et al. (1997) and in 
the X-ray from Laor et al. (1997). The UV-optical spectrum has been scaled 
to match the flux of 3C 273 at 1000 A, near the break found in the composite 
spectrum (1050 A rest) and in the 3C 273 spectrum (919 A rest). The X-ray 
spectrum has been scaled to the UV-optical composite using the mean a ox for 
the Laor et al. sample. The dotted and dashed curves represent the spectral fits 
to 14 observations of 3C 273 over a 30-day period with ROSAT (LMP95). All 
but two of these observations (dashed) gave a m p = 1.64—1.78. 

Finally we have also plotted the results of an extensive set 
of observations of 3C 273 with ROSAT in 1992 December to 
1993 January by LMP95. These authors obtained spectra at 2- 
day intervals over a period of a month. They found that the 
only successful fits to their spectra were obtained with a sum 
of two power laws, where they constrained the hard compo- 
nent to have a/, = 0.5 as determined by EXOSAT and Ginga 
2-10 keV observations, and also constrained the interstellar ab- 
sorption at low energies to be given by the Galactic value of 
Nh = 1.84 x 10 20 cm" 2 . Their soft component was then found 
to have a s = 1.7 for all but two of their observations, where 
they found a s = 1 .4. We have plotted the fits given in Table 6 of 
LMP95 as dotted-line curves in Fig. 6, with the two discrepant 
observations as dashed-line curves. It is immediately apparent 
that these observations of 3C 273 match quite well the compos- 
ite radio-loud X-ray spectrum of Laor et al. (1997) as we have 
scaled it to match the UV-optical composite, which itself was 
scaled to match our observation of 3C 273 at 1000 A. They are 



also seen to extrapolate very well to match the extreme ultra- 
viole t part of our composite UV-optical spectrum (Z heng et al. 
1997). 

The upper panel of Fig. 6 thus demonstrates the ultraviolet 
peak in the energy distribution of a single quasar, 3C 273, and a 
likely power-law connection extending from the Lyman limit to 
the soft X-ray region, with otEuv = 1 .7 ± 0.36, though a substan- 
tial extreme ultraviolet gap necessarily remains unfilled. The 
lower panel of Fig. 6, on the other hand, makes use of compos- 
ite quasar spectra constructed in the UV-optical and the X-ray 
regions (suitably scaled) along with soft X-ray data for 3C 273 
(LMP95), to reveal a very similar spectrum, with a significantly 
smaller extreme ultraviolet gap remaining unfilled. 

The striking similarity between the 3C 273 spectrum and the 
spectral composite provides genuine physical support for the 
reality of the composite spectral shape derived by Zheng et al. 
(1997). (We note that 3C 273 made no contribution to the shape 
of the composite as only quasars with z > 0.33 were included 
in the sample used.) One can argue that composite spectra in 
and of themselves are unphysical. Zheng et al. (1997) note that 
there is a wide dispersion in spectral indices among individual 
quasars. Koratkar & Blaes (1999) show several examples of 
widely differing QSO spectral shapes. So, there is no guaran- 
tee that the composite assembled by Zheng et al. (1997) actu- 
ally represents a real spectrum. The fact that the spectrum of a 
single object such as 3C 273 exhibits such remarkably similar 
characteristics bolsters the case for the broad applicability of 
the composite spectrum. Taken together, the composite spec- 
trum and the 3C 273 spectrum provide compelling evidence for 
a power-law spectrum from the Lyman limit to the soft X-ray 
band in quasars, with a typical value of a = 1.7-2.2. 

5. ACCRETION DISK MODELS 

To place our empirical results on the continuum spectral 
shape of 3C 273 in a physical context, we compare our spectra 
to accretion disk models. Such models are well suited to pro- 
ducing spectra that peak in the far-UV, and the spectral break 
near the Lyman limit in the quasar rest frame suggests a mech- 
anism related to the large change in opacity at the intrinsic 
Lyman edge. The presence of a break rather than an edge, 
however, plus the seeming extrapolation of the high-frequency 
power law to soft X-ray energies are suggestive of the appear- 
ance of an accretion disk spectrum with an intrinsic Lyman 
edge feature that has been Comptonized by an external medium. 
Comptonization (and relativistic effects in the disk) smear out 
any intrinsic features in the disk spectrum as well as add ing a 
power-law high-energy tail (Czerny & Zbyszewska 1991). As 
we will note below, however, the observed break is still sharper 
than what can be accommodated by our simple, semi-empirical 
models. 

To fit the broad-band spectral shape, we compute disk spec- 
tra in the Schwarzschild metric that are a sum of blackbody 
spectra representing rings in the disk running from an inner ra- 
dius of 6 gravitational radii (fc = GMbh/c 2 ) to 1000 tq. The 
blackbody spectra are modified by an empirical Lyman-limit 
feature as described by Lee, Kriss, & Davidsen (1992) and Lee 
(1995). For the Comptonization we assume that a hot spher- 
ical medium with an optical depth to Compton scattering of 
r e = 1 surrounds the disk. We use the formulation of Czerny 
& Zbyszewska (1991) as described by Lee, Kriss, & David- 
sen (1992) to calculate the effects of the Comptonization. The 
spectral index of the high energy tail in this calculation depends 
almost entirely on the Compton y parameter, y oc r 2 J c , where 
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T c is the temperature of the Comptonizing medium ( Lee 1995 ), 
so the free parameters in our accretion disk model are the mass 
accretion rate, m, the mass of the central black hole, Mbh, the 
optical depth at the Lyman edge, r Ly , the inclination of the disk, 
i, and T c (since we keep r e fixed at 1). 

Our fits use the same wavelength intervals as the power-law 
and broken power-law fits in §3.2 and the same emission and 
absorption components. All fits have the extinction fixed at 
E(B-V) = 0.032. Since there is little data in our HUT spectrum 
shortward of the spectral break at ^920 A in the rest frame, it 
is the slope and intensity of the soft X-ray spectrum that largely 
determines T c . This has a best fit value of T c = 4.1 x 10 8 K in 
all our fits, regardless of the other parameters. The shape of the 
intrinsic disk spectrum is invariant for a fixed ratio of rhjM\ H \ 
this ratio is determined largely by the wavelength of the peak 
in the spectral energy distribution. The normalization is deter- 
mined by the mass accretion rate m. For conversion of flux to 
luminosity, we assume a Hubble constant Hq = 75 km s" 1 Mpc" 1 
and a deceleration parameter qo = 0. The strength of the spectral 
break in our model is related both to the peak in the spectral en- 
ergy distribution as well as the optical depth at the Lyman edge, 
Ti y . The break wavelength and its sharpness constrain the best 
fit inclination i. (Higher inclinations lead to breaks at shorter 
wavelengths that are more smeared out.) 

Our best fit to the Astro- 1 spectrum, for which simultaneous 
X-ray data and nearly simultaneous near-UV and optical data 
are also available, yields \" = 1609 for 1622 data points and 
82 free parameters. The best-fit parameters are summarized in 
the first column of Table 8. (The emission and absorption line 
parameters have nearly identical values to those presented from 
our previous empirical fits.) The quoted errors are 90% confi- 
dence for a single interesting parameter, and they correspond to 
A\ 2 =4.6. The Comptonized accretion disk fits significantly 
better than a single power law at the nominal value for Galac- 
tic extinction, but worse than the best fitting broken power law. 
The sharp break present in the spectrum is difficult to match 
with the smoother shape of our simple accretion disk model, 
even with a Lyman-limit feature. 

TABLE 8 

Accretion Disk Models for 3C 273 



Parameter 



Astro- 1 



Astro-2 



M BH (10 s M ) 


7.1 ±0.3 


12 ±0.4 


m(M Q yr 1 ) 


13±0.13 


12±0.12 


i (deg) 


60_ 10 




Thy 


0.5 ±0.5 




Ts (10 8 K) 


4.1±0.2 


4.1±0.2 


E(B-V) 


0.032 


0.032 


X 2 /dof 


1609/1540 


1752/1497 



The best fit to the Astro-2 spectrum is also summarized in 
Table 8. The most significant difference between the Astro- 1 
and the Astro-2 fits is that since the Astro-2 spectrum is redder 
overall, its fit favors a higher mass black hole. In all other re- 
spects, however, the Astro- 1 and Astro-2 disk fits are similar. 
Reflecting the nearly identical flux levels of the two observa- 
tions, the mass accretion rates are comparable. Both fits favor 
an inclination i = 60°, and the Lyman-limit optical depth in each 
is low. For Astro- 1, Thy = 0.5, and Thy = 0.0 also gives an accept- 
able fit; for Astro-2, Thy = 0.0. Thus in the context of our simple 
accretion disk model a Lyman edge feature is not formally re- 
quired to achieve an acceptable fit to our data for either Astro- 1 
or Astro-2. However, the fact that a broken power-law actually 



provides a better fit, and that the break is very close to the Ly- 
man limit, hints that there probably is an effect associated with 
an opacity change at the Lyman limit. We suggest that future 
more sophisticated disk models should endeavor to match this 
feature. 

Fig. 7 shows the best-fit accretion disk spectrum in compar- 
ison to the HUT Astro-1 data and longer-wavelength points. 
Here one can see that the rather smooth disk model does not 
fully accommodate the sharp spectral break at ^1000 A. In 
the near-UV range the accretion disk continuum falls below the 
data largely due to Balmer continuum and Fe II line emission 
that is not included in our model. At longer wavelengths, the 
accretion disk falls below the data even more as we have not 
included any emission components (either a power law or dust 
emission) that are mainly responsible for the near-IR contin- 
uum. Fig. 8 again shows the best fit over the whole frequency 
range from the hard X-ray to 7000 A in the visible. The Comp- 
tonized tail of the accretion disk spectrum provides an excellent 
match to the soft X-ray excess present in the ROSAT and the 
BBXRT data. 
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FIG. 7. — Quasi-simultaneous UV-optical spectrum of 3C 273 from 1990 De- 
cember, including data from HUT, IUE, and the KPNO 2. 1 m telescope (kindly 
provided by R. Green), corrected for reddening with E(B— V) = 0.032. The data 
have been smoothed with a 9-pixel boxcar filter. The smooth curve is the best- 
fitting accretion disk model as described in Section 5, with parameters given in 
Table 8. 

Several aspects of our accretion disk fits illustrate the short- 
comings of our simple model. Our best-fit inclinations of 
i = 60° are a bit high given the superluminal jet in 3C 273 that 
suggests the disk normal should be directed close to our line 
of sight. However, inclinations in superluminal sources can be 
as high as 30-45°, and this provides an acceptable fit for our 
Astro-2 observations. An inclination of 0° is excluded at high 
confidence in both data sets (Ax 2 = 27 for the Astro-1 data). 
Also, as mentioned earlier, the inclination is determined mostly 
by the wavelength of the spectral break. A Kerr disk would 
produce bluer breaks that are more smeared out at much lower 
inclinati on than in a Schwarzschild disk (e.g., see L aor & Net- 
zer 1989). 

Both the Astro-1 and Astro-2 fits have accretion rates that 



10 



FAR-UV PEAK OF 3C 273 



m 
I 



o 

00 

u 

CD 



o 



-9 



-10 - 



-11 



T 



~ i 1 1 r 



T 



IUE 



KPNO 



HUT 



3C 273 
December 1990 




15 



16 17 

l0 g ^rest ( HZ ) 



18 



FIG. 8. — Quasi-simultaneous broad-band spectrum of 3C 273 as described in Fig. 5, with the best-fitting Comptonized accretion disk model added to a hard 
X-ray power-law with a = 0.5. The sharp turn downward of the observed data below the model at the shortest ultraviolet wavelengths is simply absorption by the 
converging hydrogen Lyman series, and it is fully accounted for in our models (see Fig. 3). 



exceed the thin-disk limit of L/L£dd ~ 0.3 for a Schwarzschild 
metric. For the Astro-1 fit, L/LEdd = 0.46; for Astro-2, 
L/Lgdd = 0.25. A Kerr metric would also alleviate this prob- 
lem: the limit for a thin disk in this metric is higher, and, at 
the lower inclinations it would require, the accretion rate would 
also be lower. 

Finally, while the Astro-1 and Astro-2 fits have similar mass 
accretion rates, their best-fit black hole masses differ by nearly 
a factor of two. This illustrates a fundamental difficulty of 
simple steady-state accretion disk models in dealing with vari- 
ability: it occurs on time scales more rapid than the viscous 
time scale that governs the applicability of steady-state mod- 
els. In a steady-state model for a given object, the only vari- 
able that can vary freely is the mass accretion rate. There- 
fore, one should observe correlated changes in flux and color, 
with brighter states exhibiting bluer disk spectra. However, 
the Astro-1 and Astro-2 data are nearly identical in flux, but 
they have significantly different colors. Without a good model 
for accommodating variability in disk models, one can see that 
factors of several uncertainty in the actual physical parameters 
are easily present in our results. It is reassuring to note, how- 
ever, that our values for the black hole mass based on accretion 
disk models (7.1 - 12 x 10 8 M Q ) bracket the independent value 
of 7.4 x 10 8 M Q obtained by Laor (1998) using reverberation- 



mapping data. 

We conclude that our accretion disk models of the Astro-1 
and Astro-2 spectra provide an overall qualitative physical char- 
acterization of the 3C 273 spectrum. While the models have 
quantitative shortcomings, they provide an empirical guide for 
potentially more sophisticated models. The most important as- 
pects of our simple description in terms of blackbody emission 
and Comptonization are that it (1) accounts for the peak of the 
ultraviolet spectrum, (2) produces a spectral break in the far- 
UV, and (3) has a hard power-law tail that extends to the soft- 
X-ray band and matches the observed soft-X-ray spectrum. 

6. SUMMARY 

We have presented absolutely-calibrated ultraviolet spec- 
trophotometry over the 900-1800 A range for the quasar 
3C 273, obtained with the Hopkins Ultraviolet Telescope on the 
Astro-1 mission in December 1990 and on the Astro-2 mission 
in March 1995. In both observations the continuum displays a 
change of slope near the Lyman limit in the quasar rest frame. 
At longer UV wavelengths the continuum is well-represented 
by a power-law of ct\ = 0.5. Shortward of the Lyman limit, 
however, the continuum slope has a-i = 1.7 ±0.36, where the 
uncertainty includes our uncertainty about E(B-V) at the level 
of ±0.01. 
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The energy distribution per logarithmic frequency interval 
vf v therefore has a peak close to the quasar Lyman limit. The 
short wavelength extreme UV power-law extrapolates very well 
to match the soft x-ray spectrum of 3C 273 obtained nearly si- 
multaneously in the case of Astro- 1 with BBXRT and ROSAT. 
The soft x-ray data themselves give a s = 1.7(±0.1) (LMP95), 
so the combined UV and X-ray data suggest auv-x = 1-7. 
While some models for the photoionizing extreme ultraviolet 
radiation from quasars have a peak — the so-called "big blue 
bump" — at a wavelength much shorter than the Lyman limit 
at 912 A with a peak flux much higher than what is actually 
observed at both lo nger and shorter wavelengths (e.g., M ath- 
ews & Ferland 1987; IBechtold et al. 1987t Gondhalekar et al. 
1992), we find that in 3C 273 the peak of the big blue bump 
occurs very close to the Lyman limit. The photoionizing flux 
from 3C 273 is given by f v = fuSy jvi£f xl extending from the 
Lyman limit to 1 keV, although there is still a gap in the data 
covering one decade of frequency in the extreme ultraviolet. 

Analys is of composite qu asar spectra con structed in the ul - 
traviolet flZheng et al. 1997| ) and the X-ray ( |^aor et al. 1997| ) 
bands leads to the same result. The similarity between the spec- 
trum of a single object such as 3C 273 and these composites 
lends physical credence to their applicability. It therefore ap- 
pears to be true of quasars in general that the peak of the con- 
tinuum energy distribution occurs near the Lyman limit, and 
that the ionizing continuum is well-represented by a power-law 
of energy index auv-x = 1.7-2.2 extending to about 1 keV, 
where a separate hard X-ray component begins to dominate for 
higher energies. The "extreme ultraviolet gap" in the composite 
quasar spectrum is now only half a decade in frequency, how- 
ever. This can be reduced somewhat further by extending the 
work of Zheng et al. (1997) to include more observations of 
high redshift quasars in the far ultraviolet band. 

The general shape of the optical, ultraviolet, and soft x-ray 
spectrum of 3C 273 is consistent with that of an optically thick 
accretion disk around a massive black hole, which is itself sur- 
rounded by a hot medium that modifies the emergent spectrum 
through the inverse Compton effect, as suggested by Czerny & 



Zbyszewska (1991), Lee et al. (1992) and Lee (1995). In our 
model the spectral break is due to the thermal peak of the accre- 
tion disk spectral energy distribution with a minor contribution 
from Lyman edge absorption. Both features are broadened and 
blurred by Comptonization in a surrounding hot medium, per- 
haps a corona or wind emanating from the disk. Comptoniza- 
tion of the thermal photons from the disk also produces the 
power-law ta il extending to the soft X-ray region (e.g., C zerny 
& Elvis 1987 ; [Maraschi & Molendi 199o[ R oss, Fabian, & Mi- 
neshige 1992). Our simple model is unable to fully account for 
the sharpness of the spectral break, however, and we suggest 
that a more sophisticated treatment of the Lyman-edge region 
may eventually provide a better match to these data. 

Assuming a Schwarzschild black hole with an inclination 
of 60 degrees, the UV spectrum is fit with a black hole mass 
Mbh = 7 x 10 s Mq with an accretion rate M = 13 M Q /yr. Su- 
perposed on the disk spectrum is an empirically determined Ly- 
man edge with Ti y = 0.5. The hot corona or wind is required to 
have a Compton parameter y rj 1, which is obtained, for ex- 
ample, with T es = 1 and T e = 4 x 10 8 K. A temperature of this 
order appears plausible for a corona where Compton cooling 
balances Compton heating by the hard x-ray flux in 3C 273. 
These results are quite similar to those we found previously by 
fitting simi lar models to a composite quasar spectrum (Z heng 
et al. 1997). There we found, assuming a typical inclination 
of 30 degrees, a black hole mass Mbh = 1.4 x 10 9 M Q and an 
accretion rate M = 2.8 M & /yr. While the best-fit parameters 
would change somewhat in the more realistic case of a rotating 
(Kerr) black hole (which we have not calculated), we believe 
these results provide strong evidence in favor of the hypothesis 
that quasars are powered by accretion onto massive black holes. 
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